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Abstract 

We generalize the worldline variational approach to field theory by introducing a trial action 
which allows for anisotropic terms to be induced by external 4-momenta of Green’s functions. By 
solving the ensuing variational equations numerically we demonstrate that within the (quenched) 
scalar Wick-Cutkosky model considerable improvement can be achieved over results obtained pre¬ 
viously with isotropic actions. In particular, the critical coupling associated with the instability of 
the model is lowered, in accordance with expectations from Baym’s proof of the instability in the 
unquenched theory. The physical picture associated with a different quantum mechanical motion 
of the dressed particle along and perpendicular to its classical momentum is discussed. Indeed, 
we find that for large couplings the dressed particle is strongly distorted in the direction of its 
four-momentum. In addition, we obtain an exact relation between the renormalized coupling of 
the theory and the propagator. Along the way we introduce new and efficient methods to evaluate 
the averages needed in the variational approach and apply them to the calculation of the 2-point 
function. 


1 Introduction 


Feynman’s variational method for evaluating functional integrals, well known in condensed matter 
physics 1]^, makes use of a trial action containing variational parameters. Even though the exact 
action of the theory under consideration can in general be quite complicated, the trial action itself 
can be at most quadratic in the degrees of freedom of the theory because the only functional integrals 
which can be performed analytically are Gaussian ones. Clearly this imposes a significant constraint 
so that one might at first sight suspect that such a variational method, even if valid in principle, is 
likely to be of limited use in practice. This turns out not to be the case, however: Feynman formulated 
the variational calculation on the level of a non-local action for the polaron, obtained after integrating 
out the phonon degrees of freedom. This non-locality, which is shared by the trial action used by 
him, introduces considerable freedom for the variational principle to work with and hence he obtained 
numerical results for, e.g. the polaron’s groundstate energy, which differed ~ for a very large range of 
coupling constants - from the exact ones by at most 2 % percent. 

Encouraged by this success we have extended this method to relativistic field theory [Q| - fol¬ 
lowing the pioneering, but largely forgotten, work by K. Mano Q. In these investigations Eeynman’s 
approach was applied to a relativistic scalar field theory of two fields ip and <h with Lagrangian (in 
Minkowski space) 

-C = ^ + i + (1-1) 


where we refer to the light field ip as “meson” and to the heavy one (<h) as “nucleon”. This theory 
is the Wick-Cutkosky model |9|, 10| and is of interest as it models, in a rather rough fashion, a 
simplified nucleon-pion theory without the complications of spin and isospin (and chiral symmetry). 
Due to the close similarity to the polaron model (indeed, the dressed nucleon might be referred to 
as “relativistic polaron”, the only essential difference being the space-time dimensionality) one would 
expect the variational method to work equally well in this setting. Despite of this similarity, there 
are noticable qualitative differences: most importantly, while the polaron is a stable (quasi-) particle 
at all couplings, it turns out (see Refs, ^]) that the variational equations of the Wick-Cutkosky 
model only have (real) solutions below a critical coupling of ac ~ 0.815, where 


a = 


g 


dvrM^ 


( 1 . 2 ) 


is the dimensionless coupling constant and M the physical mass associated with the field d*. The 
existence of a critical coupling is likely to be a remnant of the well-known instability of cubic scalar 


field theories |11]. 


There has been renewed interest in this theory and its instability in recent years. Ahlig and 
Alkofer [12| have analysed the stability of the model using Dyson-Schwinger equation methods and 
have found a critical coupling rather similar to the one mentioned above. Tjon and coworkers 
have studied the theory (without the radiative corrections which are responsible for the instability) 
using Monte-Carlo techniques and have recently also alluded to evidence of a critical coupling |^. It 
has also been suggested by these authors [IT] that the apparent stability of the quenched theory below 
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some critical coupling is a real effect []and not just a byproduct of approximations [^]. In principle, 
the “true” ac in the quenched theory could be found numerically by Monte-Carlo simulation similar 
to the way the true ground state energy of the polaron has been determined by stochastic methods 
[ll^ . Notwithstanding the interest in this critical coupling, we also note that there has been a great 
deal of interest in recent years in field theories with imaginary or negative coupling constants (which 
would lead to stability even for the unquenched theory) which, although they have non-hermitian 
Hamiltonians, are nevertheless argued to have a positive-definite spectrum by virtue of their PT- 
symmetry []^|. Finally, of course, there is continued interest in using this field theory as a vehicle for 
examining boundstate problems []l^ . 

In recent applications we have treated a more realistic fermionic theory, viz. Quantum Electro¬ 
dynamics, by worldline variational methods and obtained a compact non-perturbative expression 
for the anomalous mass dimension. This shows that this approach successfully describes the short- 
distance singularities of relativistic quantum field theories. It is remarkable that at the same time 
the long-distance behaviour is also caught to a large extent as evidenced by the proper threshold 
behaviour of scattering amplitudes in the scalar model S|] or the correct exponentiation of infrared 
singularities in QED 1^ . 

Encouraging as these specific results are, one would like to assess the reliability of the variational 
results in general. There are two obvious ways one can do this: the approximation at the heart of 
Feynman’s variational calculation is the cumulant expansion and corrections can be systematically 
calculated, as was done for the polaron in Ref. Alternatively, one can make use of more general 

trial actions, thus allowing the variational principle more freedom to work with. It is the latter course 
of action which we pursue here, concentrating for simplicity on the simple scalar field theory defined 
in Eq. O) and considered in Refs. |Q, p. The results presented there employed a number of different 
forms for the trial action, culminating in the most general (isotropic), non-local, quadratic trial action 
possible. It might seem, because the trial action is required to be quadratic, that it is hard to improve 
on this. However, as has been explored within the context of the polaron in Ref. |^|, one can make 
use of the external direction provided by the particle’s four-momentum to construct more general, 
anisotropic, trial actions. Eor the polaron this extra freedom only yields marginally improved results 
due to the nonrelativistic nature of that problem |^]. One might expect however, and we shall indeed 
find, that significant improvements can be achieved in the present, relativistic, theory. 

The organization of this paper is as follows. In the next Section, after briefly summarizing the 
essential details of the variational method (for more details, we refer the reader to Ref. |Q), we describe 
the anisotropic trial action and derive the relevant variational equations. In Sec. ^we present numerical 
results and in Sec. 0 we conclude. Technical details are collected in five appendices. 


^In our opinion, however, this suggestion cannot be correct as it conflicts with standard expectations based on the 
behaviour of high orders of the perturbative expansion of the theory: The contribution from the order to the 
Euclidean coordinate space 2-point function grows roughly like the number of Feynman diagrams at that order. This 
growth is factorial in both the unquenched and quenched theory. Moreover, as the expansion is in rather than 
g, the contribution from each order is of the same sign. Therefore the series is not Borel summable, indicating that a 
perturbative expansion is an expansion around the wrong vacuum. Quantitative calculations summarized in Appendix^ 
confirm this expectation. Also, it seems that the “proof” in is flawed inter alia by wrongly asserting that the true 
mass is larger than the variationally estimated value which invalidates all subsequent steps. 
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2 The variational method 

2.1 Feynman’s variational method applied to the Wick-Cntkosky model and the 
anisotropic trial action 

In the polaron variational approach the key elements are 


i) The reduction of degrees of freedom by using the particle (“worldline”) representation and by 
integrating out the mesons/photons. For example, in the one-nucleon sector of the scalar theory 
defined by Eq. O) this is achieved by using the Feynman-Schwinger representation (as well as 
neglecting the functional determinant if one is working in the quenched approximation) 


1 

—— Mq + 2g(p{x) + iO 


i 

2ko 


dT exp 


i 

2ko 



^0 + T 


( 2 . 1 ) 


Here kq is a free (positive) parameter which reparametrizes the proper time T without affecting 
the physics The exponential factor on the r.h.s. of Eq. ( p.lj ) can be considered as a proper¬ 
time evolution operator and may be represented as a quantum mechanical path integral. With 
(^(x) now being a c-number quantity the integration over the mesonic field can be performed 
resulting in an effective action for the four-dimensional trajectory x{t) of the nucleon 


S[x] 



^ L 


pto-\-T 


I to 


dt2 


d^k 


(27r)^ -|- zO 

• exp {-i/c • [x(ti) - x(t 2 )] } = 5o -I- S'! . (2.2) 


Here the starting time to is free since only the proper time interval T matters. Convenient 
choices are to = 0 or to = —T/2. 

ii) The use of the Eeynman-Jensen variational principle 

j Vxe"-^^ exp[i{S - St)] = Vxe^^^'^ ■ {exp[i{S - St)])st 

~ T>xe*'^‘^ • exp [i(S'- , (2.3) 

where < ... >St refers to averaging with weight function (normalized to < 1 >St— 1) 

~ indicates equality at the stationary point of the r. h. s. under unrestricted variations of the 
trial action St- 


iii) A quadratic two-time trial action 


Sr (1) = 5o + 


^ k 


dtidt 2 fih - t 2 ) x{ti) - x{t 2 ) 


l2 


(2.4) 


in which the non-quadratic terms in the true action S are approximated by an even retardation 
function f{cr). Since to is arbitrary (proper time-translation invariance) this function can only 

^It also allows an easy (formal) switch to Euclidean space by setting kq = ike and changing the sign of 4-vector 
products. Proper times remain unchanged and the action S becomes iSE- 
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depend on the time difference a = ti — t 2 . Note that the true action as well as the trial action 
are invariant under a constant translation in x. Therefore one also could take 


-^iso (2) _ 


rtQ-\-T 


= So- zkq / dtidt 2 gih - t 2 ) x{ti) ■ x{t 2 ) , 
Jto 


(2.5) 


because two integration by parts and neglect of boundary terms show that this is equivalent to 
Eq. (U) with /(cr) = g{a). 

Being mainly interested in the on-shell-limit of the Green functions of the theory it is convenient to 
include the Fourier transform over the endpoints into the path integral (“momentum averaging” in 
the parlance of Ref. [|| ) and to consider 


5[x] = p - X + 5[x] 


( 2 . 6 ) 


The nucleon’s propagator may be calculated from the generating functional and the worldline repre¬ 
sentation via 


G2ip) = 


1 


2iKo 


dT exp 




/ Vx 
f Vx 


(2.7) 


-■92{p,T) 

Note that we have normalized the functional measure by dividing by the path integral involving the 
free action So (i.e. Eq. (|2.6| ) with 5 = 0). As mentioned, the path integrals in Eq. (P) include an 
integral over the endpoint x - hence the tilde over the x. 

The variational approximation to this propagator, making use of Jensen’s stationary principle, has 


the form (2.7) but with 


52 (p,r) gTiP^T) = 


fVx 
f Vx 


exp 


. / Vx (5[x] — St[x]) 
f Vx 


( 2 . 8 ) 


Of course, when the extended trial action St equals S this expression gives the exact propagator in 
Eq. (2.7). On the other hand, as long as St is quadratic in x{t), all functional integrals appearing in 
Eq. (|2.8D may be performed analytically. The most general linear -|- quadratic isotropic trial action, 
respecting symmetries such as translational invariance, may be written as 


srix] = xp-x+sr 


(2.9) 


where A is an additional variational parameter. Note that the structure of the linear term is basically 
fixed by requiring that the action is a scalar and by time translation invariance: there is only one ad¬ 
ditional four-vector, the external momentum p, with which the velocity four-vector can be contracted 
and any modification of the free term p ■ x = p ■ f dtv{t) can only be done by multiplication with a 
constant, but not with a function of the proper time. In contrast, the retardation functions / or 5 in 
the quadratic term must be functions of the time difference a = ti — t 2 - 

The key to the success of Feynman’s variational method is that the different functional dependence 
on the path x{t) in the exact action in Eq. (2.6) as compared to the trial action in Eq. ( p.4| ) can be 
compensated for by the variational retardation function f{cr). For example, in the path integral in 
Eq. (13) contributions from [x(ti) — x(t 2 )]'^ ~ 0 are greatly enhanced because of the UV divergence 
in the interacting part of the action Eq. ( |2.6D . In the path integral in Eq. (^) the equivalent 
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enhancement is provided by a divergence of /(cr) as cr —> 0. Similarly, infrared physics is simulated by 
the large-cT behaviour of this function. However, in general it is to be expected that in the functional 
integral over exp(iS'[x]) paths which contain segments — x^(t 2 ) which are parallel to the 

momentum will receive a different weighting to those containing segments which are perpendicular 
to p^. Because f{a) is scalar, the trial actions (^, cannot differentiate between these two. 
Rather, the behaviour of the retardation function /{cr) will encapsulate some compromise of this 
directional information. 

The situation is different, however, if we generalize the trial action by explicitly making use of the 
external vector The most general covariant, quadratic, anisotropic trial action may be written as 


^aniso 1 





pP‘pl'\ 






Xy{tl) 


Xu{t2) 


( 2 . 10 ) 


In this action we now have two independent retardation functions fT{<x) and /l((t) and we shall show, 
in the following, that these indeed encode just that directional information in the functional integral 


which we have discussed above. Equivalently, just as in the isotropic case (i.e. Eq. (|2.5D ), one could 
take 


caniso2 _ q .*2 

bt — *bo - 


^to+T 


ho 


dtidt2 


I , ! \ ( ay P^P’' 

gUcr) —^+ 5 t ( o -) 5 ^-^ 


X^j,{tl)x„{t2) ( 2 . 11 ) 


instead. 


2.2 Mano’s equation and functional averages 


In Appendix ^ it is detailed how the averages required in Eq. (^^) may be calculated easily by a 
method which is more transparent and efficient than the Fourier series expansion used in Ref. |0. In 


particular, Eq. (B.2C) demonstrates that all trial actions may be written as 




( 2 . 12 ) 


if expressed in terms of velocities v. Allowing for a Lorentz structure of A, i.e. (u|A|u) = 
this is obviously the most general linear + quadratic trial action. In principle one could work out the 
required averages with this trial action and let the variational principle determine the form of A and 
b. That this indeed leads to an anisotropic trial action is sketched in Appendix ^ However, since 
there is only one four-vector available for the 2-point function, namely the external momentum p, it is 
clear that b must be proportional to p and that A must have the decomposition ( B.19| ), which reflects 
the Lorentz structure of the retardation functions in Eqs. (T 1, iH). 

In the following we will only consider the limit of the 2-point function where —> M^, i.e. the 


on-mass-shell limit. It is well known |24| that in order to obtain the pole and the residue of the 
nucleon propagator 

p2 - M2 + iO ’ 

the proper time T has to tend to infinity. To obtain the physical mass in terms of the bare mass only 
the leading large-T term in the exponential of §2 is of relevance (see Eq. (D.lt) of Appendix 
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Performing the T-integration in Eq. (2.7) and setting gives 


Ml = {2X-X^)M^ -2{n + V) 


(2.14) 


This equation, whose content is discussed in detail below, was termed ‘Mano’s equation’ in Ref. |^. 
There it was also shown that for fixed physical mass M (and fixed UV regulator) the variational 
principle in fact provides a lower bound on the exact bare mass, i.e. 


MQ(exact) > Mq (Eq. (2.14)) 


(2.15) 


The parameter A in Mano’s equation (2.14) is related to the original variational parameter A by 

A 


A = 


AIl(O) 


(2.16) 


As shown in Eq. (p.8 ) those averages not involving the coupling a are contained in the ‘kinetic term’ 
XI, which in d dimensions is given by 


Vl[At,Al\ = H- 2 —^[^r] 


with 




dE 


log Al,t{E) + 


1 


- 1 


(2.17) 


(2.18) 


Al,t{E) 

Note that when At = At this reduces to the isotropic result and that the factor d — 1 in front of the 
first term arises because in d dimensions there are d — 1 dimensions transverse to p^. Erom Eqs. ( D.IO , 
D.ll ) it is seen that the transverse and longitudinal profile functions At,t are essentially the cosine 
transforms of the transverse and longitudinal retardation functions, respectively. Since log A + 1/A — 1 
is positive for A > 0 the kinetic term provides the restoring force to balance the attractive interaction. 

The ‘potential’ term, essentially the average of the interaction dependent part of the action, is 
given by Eq. (D.15) which in d = 4 dimensions reads 


V = - 


Svr^ 


da / du 

Jo 




-exp 

a 


exp 


2ko 


m ^i(cr) 


1 — u 


* 2 
- m a - 

2kq u 


1 — u 
u 

f_ 

87r2 


+ 


X'^M^a^ 




-u 


- — 


(2.19) 


In addition, we have used a Pauli-Villars regulator as in Ref. to make the proper time integral 
convergent for small a. As in Refs. |Q, ^ the divergent logarithm on the right hand side will be taken 
over to the left hand side of Eq. (2.14), with the replacement of Mq by the finite 


A 2 

Ml = Ml-j-^\og\ 

47r^ 




( 2 . 20 ) 


The essential new element in Eq. ( 2.19|) as compared to the isotropic case is the appearance of the 
anisotropy factor 

Vr(o-) 


FLT{u,a) = 


1 + 




- 1 U 


-1 


( 2 . 21 ) 
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Also, note that there is now one ‘pseudotime’ /r^((T) for each of 




dE sin^ ^ 


vr Jo Al,t{E) 


( 2 . 22 ) 


This function was named as such in Ref. @ because, as can be seen rather straightforwardly from 
Eq. ( 2. 221) , it is proportional to cr (fi > 0) for small a as well as for large a (and, indeed, is equal 
to a when the interaction is turned off, i.e. when A{E) = 1). Actually, these pseudotimes also have 
a direct physical interpretation: Let us consider, initially, the average separation x{ti) — x{t 2 ) when 
weighted with the exponential of the trial action. Using Eqs. (B.30, D.12| ) this is easily calculated 
and one obtains 


\p^ 


o 


(2.23) 


(xnti)-x^(t2))5, = - 

The meaning of this expression is clear: the average (four-) displacement that the nucleon undertakes 
in the proper time interval a = ti — t 2 is given by its (four-) velocity • a. It is natural then to view 
kq/\ as an effective mass |^. That this interpretation is not at all unreasonable was seen for the 


3-dimensional polaron problem discussed in Ref. |23|, where the variational equations for A directly 
lead to this result. 


It is very instructive to consider now the mean square displacement. With Eq. ( D.13 ) one obtains 


{[X^{E)-X^{t2)] K(tl)-x"(t2)])5* = 


U U 

rfliOi - 


Kq 


- 0.2 


aVp" 




KO 


. (2.24) 



Figure 1: The physical interpretation of the pseudotimes The dotted circle denotes the 

noninteracting particle with quantum mechanical uncertainty in (x^) given by Brownian motion alone. 

The first term is again the contribution from the straight line (i.e. classical) motion. The last 
term characterizes both the quantum mechanical deviations from this as well as the ‘jiggling’ the 

^ Eq. ( ^.2[ ) shows that kq is the “mass” of the particle in the worldline description. Integration over the proper time 
T with the weight exp(—iMor/(2Ko)) gives the particle its actual (bare) mass. Note that in the nonrelativistic limit 
the proper time can be identihed with the ordinary time when kq = M is chosen |p5||. 
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nucleon experiences because of the constant emission and re-absorption of the pions. For example, 
as mentioned above the pseudotimes reduce to a and A = 1 when the interactions are turned off. In 
this case the first term in Eq. ( |2.24 ) represents a constant drift and the second just corresponds to 
the result well known from Brownian motion that the mean square distance (in each direction) grows 
linearly with Euclidean time (see footnote 2). Turning on the interactions modihes this Brownian 
motion, and since 

. V .... \ 

(2.25) 


f 2/ 2/ , 2( av P^P’'\ 

[P ) = Pl{^) + pA^) (^ 5 ^ - j 


p. ^ p. 

this modification may be different in the longitudinal and transverse directions. 

The results are summarized in the cartoon in Eig. 2, anticipating the numerical results from 
Section ^ Note that the possibility of distortion in the longitudinal direction is only present because 
the trial action has been allowed to be anisotropic. It is shown in Appendix that for nonvanishing 
momentum p the variational principle, if given sufficient freedom, demands such a structure. 


2.3 The variational equations 


By demanding that the independent variation of Mano’s equation (ie. Eq. (p.l4 )) with respect to 
the parameter A and the profile functions A^iE) and Ai[E) must vanish, we obtain the variational 
equations for these functions. Since a numerical solution is only feasible in Euclidean space we choose 
kq = Ike- After some work one obtains 


I ^2 I 


^— [ dcr f [ duuFli^{u,a) ■e{u,a) 
ke Jo PlA) Ao 


(2.26) 


At{E) 

5V 

dPricr) 


= 1 + 


8 


3keE'^ 

3g2 1 


da 


6V 


Ea 


sm 


167 r 2 pAA do 


5iA{a) 2 
duuE^i^{u,a) ■e(u,a) 


(2.27) 


AAE) = 




5pi{a) 


where 


1 + 

9 


8 


keE'^ Jo 

2 1 .1 


da 


5V 


Ea 


sm 


IGvr^ pAA Jo 


Splia) 

du u E^i^ {u, a 


1 - 


1 

i^E pUA 


e{u, a) = exp 


1 


m PlA. 


I — u 


+ 




e{u, a) 


(2.28) 


(2.29) 


2 ke V ' u ' plia) ^ 

Note that the compact form of the last variational equation results from an appropriate integration 


by parts in u. Neither Eq. {2.21) nor Eq. (2.28) reduce to the variational equation for the isotropic 
profile function given in Sec. IV. C of Ref. for Elt = 1 because the anisotropy factor Fit has 
been varied as well. However, both equations become identical when M = 0 and Fit = 1- This 
just corresponds to the fact that when the momentum vanishes (hence = 0) there is 

no source of anisotropy and hence no need for separate retardation functions in the longitudinal and 
transverse directions. In this case one also deduces from Eqs. (|2.28|) and ( 2.26|) that AA(0) = 1, 
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i.e. the original variational parameter A in the trial action 
interaction. 


decouples and is unaffected by the 


Finally, by comparison with Eq. (D.IC) one also sees that the retardation functions defined in 
Eq. ( 2.10| ) are essentially the variational derivatives of the ‘potential’ 




CL,T 


5V 


S^L,T 




CL = 1 , CT = 


(2.30) 


2.4 The residue 


For completeness we also evaluate the residue of the propagator in Eq. (|2.13|) . This requires the 
calculation of the next-to-leading terms in the large-T limit. As shown in Appendix p.2| only the 
interaction gives a 0(T'^)-contribution to logg 2 {jp,T) . Therefore following Sec. V. A in Ref. and 
using Eq. (p.23) we easily obtain in d = 4 Euclidean dimensions 


,7 1 

Zj = — exp 
A 


r 

Svr^ ke Jo 


da 


a 


Hiia) Jo 


du (u, a) ■ e(n, a) 


N 

T 


(2.31) 


As in Ref. |^, the variational equation ( 2.26|) for A has been used to simplify the denominator of this 
expression. It is also easily seen that Z does not depend on the value of ke- 

Apart from the obvious modifications due to the anisotropic trial action and the general reparametriza- 
tion gauge / 1, Eq. ( 2.3lJ) differs from the result obtained in Ref. |3| by the absence of the factor 
No = exp[—logA(O) -|- 1 — 1/A(0)]. As seen in Eq. (65) of Ref. HI this factor had its origin in a 
0(l/T)-correction to the kinetic term which does not occur in the present formulation where an 
even, time-translation invariant retardation function (or equivalently, an even profile function A{t—t')) 
has been employed throughout. The discrepancy is explained by recalling that - for practical reasons 
- the previous trial action was taken as a quadratic, diagonal form in Eourier space (see Eqs. (48) - 
(51) in Ref. [^). After transforming back to x-space it may be seen that this form also contains terms 
which are not invariant under time translations. The present formulation therefore not only respects 
this obvious symmetry of the action but also eliminates the spurious and awkward No-term ^ after 
absorbing A(0) in the variational parameter A, it does not appear explicitly anymore, neither in the 
pole mass nor in the residue of the propagator. 


2.5 The effective coupling 


Einally, although it is beyond the scope of this paper to discuss the vertex function in general (for 
the isotropic trial actions this was done in Ref. Q), it is rather straightforward to derive the effective 
coupling of the theory (i.e. the value of the vertex function at = 0) for an arbitrary trial action. 

We start with the (exact) worldline representation of the untruncated 3-point function, which is 
closely related to the corresponding expression for the propagator, i.e. Eq. (H3)- The only difference 
is an additional plane wave for the external pion with (outgoing) momentum and an integral 

over the proper time r at which the pion couples to the nucleon’s worldline: 


G 2 i{p,p' = P — q) = const. / dTexp 

Jo 


-Mf 

2ko 


Vx e 


—tp -X 


g dr 


/o 


JS[x\ 


(2.32) 


should be noted that all previous results for physical amplitudes are independent of Nq, as can be seen, e.g. in 
Eq. (52) of Ref. §. 
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For = 0 the integral over the time r just provides a factor of T, while the rest of the integrand is 
the same as the one for the propagator in Eq. (2.7). We therefore obtain 

d 


G 2 ,i(p,p) = 2g 


dM^ 


G2{p) 


(2.33) 


where the normalization has been fixed so that the correct free limit is obtained (see Ref. [||; the 
factor 2 arises because of the definition of the coupling in Eq. ( |1.1| ) without a factor 1/2!) and we 
have made use of the fact that the only dependence on the bare mass enters through the explicit 
exponential factor shown in Eq. ( 2.321) . One merely has to truncate the external legs off G 2 ,i{p,p') 
and multiply by the residue Z in order to obtain the vertex function. At = 0 this defines the 
effective coupling 2 f^phys • Hence the physical coupling of the theory and the propagator are related 
by the (exact) relation 


ilphys — 9 ^ 


d 


G2^{p) 


dM^ 


(2.34) 


In order to obtain the variational result for this, we write the propagator near the mass shell as 

N 


G2{p) = 


(2.35) 


M(M2) + (p2 _ m2)A ■ 

Here N is defined in Eq. ( 2.31 ) and A4(M2) _ ^2A — A^)M2 — 2(H + V) — Mq = 0 is Mano’s equation 
(B) & Note that the denominator is just the first few terms of the Taylor expansion of A4(p^) 
around = M^. The equality Ai'{p‘^)\p 2^^2 = A only makes use of the variational equation for A 
and hence is independent of whether an isotropic or anisotropic trial action is used. Substituting into 
Eq. (2.34) immediately yields the variational estimate for the effective coupling 
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ffphys 


A 


(2.36) 


As we know from previous calculations and also shall see in the next section, A is below 1 and hence 
the physical coupling is larger than the bare coupling. 


3 Numerical Results 


We have solved the variational equations ( 2.26| ) - ( 2.2j ) numerically for a variety of dimensionless 


couplings a (defined in Eq. (O)) while keeping the physical mass of the nucleon (pion) fixed at M = 
939 MeV [m = 140 MeV). As in the previous work we haven chosen the Euclidean reparametrization 
“gauge” ke = ^ which only affects the E, cr-scale of profile functions and pseudotimes, whereas 
A, ^l,t, V and therefore also the masses Mq, Mi do not depend on the reparametrization 

parameter. The numerical integrations were based on Gauss-Legendre integration, after suitable 
mapping of the infinite integration ranges to finite ones. The precision of the numerics was controlled 
by demanding that results do not change appreciably upon subdivision of the integration range. The 
equations were solved by iteration, with the convergence criterion being 1 part in 10®. In addition, as 
in the polaron case f^ , we have used the virial theorem to calculate the kinetic term in a completely 
different way: 




1 

4 ' 


= -nr + 




da a^ 


6V d (plia) 


5ni{a) da 


a 


+ 


6V d / ^j’(o') 


5p^{a) da 


a 


(3.1) 


®We need to momentarily retain Al(M^) because of its dependence on the bare mass. 
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As shown in Appendix ^ this relation relies on the fact that profile functions and pseudotimes are 
solutions of the variational equations and thereby represents a crucial test for the accuracy of the 
numerical calculation. Unfortunately, in the present case the derivative of the pseudotimes cannot be 
eliminated but has to be calculated numerically from 


dfsj^T 

da 




sin(ili(T) 

E 


{ ^ 



(3.2) 


which - due to less damping of the oscillating integrand than for the pseudotimes themselves - gives 
less stable results than for i- Nevertheless, relative agreement better than 3 • 10“^ was achieved 
for O and better than 6 • 10“® for Mi in the whole range of coupling constants. Numerical results for 
A and the values of the profile functions at E = 0, as well as the resulting values for the mass Mi and 
the residue Z, are tabulated in Table ^ for a variety of the bare couplings a. Also tabulated are the 
corresponding values for the physical coupling Ophys = aj}?. 

We draw the reader’s attention to the fact that for large a the pseudotimes fJ-'^'j'ia) behave like 
As is evident from the Table, as the coupling constant is increased the behaviour for 
these pseudotimes is rather different in the longitudinal and transverse directions, leading (for large 
a) to a strong ‘elongation’ of the dressed nucleon in the direction of (the 4-dimensional) motion and 
a ‘contraction’ in the perpendicular directions (for the largest coupling, this distortion is depicted in 
Fig-i- Clearly, the results obtained for A(0) with the best isotropic actions (see the second part of 
Table ||) were a compromise between A'r(O) and Aiit))- 



Figure 2: The prohle functions (left panels) and pseudotimes rp (right panels) for a = 0.4 


This compromise is also evident in Fig. where on the l.h.s. the transverse and longitudinal profile 
functions have been plotted as a function of E for a = 0.4. In the right hand panels the corresponding 
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Anisotropic Action 


a 

^phys 

A 

^l(O) 

At{0) 

Ml (MeV) 

Z 

0.1 

0.10566 

0.97287 

0.97688 

1.02821 

890.302 

0.96051 

0.2 

0.22472 

0.94339 

0.94996 

1.06154 

840.041 

0.91749 

0.3 

0.36159 

0.91086 

0.91781 

1.10207 

788.106 

0.86983 

0.4 

0.52357 

0.87406 

0.87796 

1.15355 

734.416 

0.81561 

0.5 

0.72456 

0.83071 

0.82554 

1.22369 

679.01 

0.75115 

0.6 

0.99842 

0.77521 

0.74794 

1.33366 

622.28 

0.66716 

0.7 

1.53567 

0.67515 

0.5675 

1.6310 

566.31 

0.5075 


Isotropic Actions 


a 

^phys 

Best Variational 

A A(0) Ml (MeV) 

Z 

^phys 

Feynman 

A Ml (MeV) 

Z 

0.1 

0.10563 

0.97297 

1.01508 

890.246 

0.96087 

0.10830 

0.96090 

890.23 

0.96090 

0.2 

0.22449 

0.94389 

1.03221 

839.785 

0.91918 

0.23663 

0.91934 

839.73 

0.91934 

0.3 

0.36051 

0.91223 

1.05202 

787.429 

0.87428 

0.39213 

0.87467 

787.29 

0.87467 

0.4 

0.51986 

0.87718 

1.07551 

732.971 

0.82521 

0.58627 

0.82600 

732.69 

0.82600 

0.5 

0.71306 

0.83738 

1.10439 

676.20 

0.77036 

0.83930 

0.77184 

675.70 

0.77184 

0.6 

0.96066 

0.79030 

1.14207 

616.98 

0.70672 

1.1923 

0.70940 

616.09 

0.70940 

0.7 

1.31472 

0.72968 

1.1972 

555.48 

0.62697 

1.7516 

0.63216 

553.93 

0.63216 

0.8 

2.06369 

0.62262 

1.3188 

493.55 

0.49284 

3.0654 

0.51086 

490.60 

0.51086 


Table 1: Top: The variational parameter A, the variational profile functions at i? = 0, the modihed 
bare mass Mi and the residue Z obtained with the anisotropic trial action. Also shown is the value of 
the effective renormalized coupling Ophys- Bottom: For comparison, the corresponding results obtained 
with the best isotropic action and with Feynman’s parametrization of the retardation function are 
collected in the second Table (see Tables III and I in Ref. Q). In all cases, the nucleon mass has been 
taken as 939 MeV and the pion mass as 140 MeV. 

pseudotimes, normalized by their free values, are shown as a function of the proper time difference a. 
Note that, curiously, the elongation in the longitudinal direction does not take place for small values 
of a. Rather, in this region of a the dressed nucleon is ‘contracted’ in all directions as compared to 
the free case. This change in sign of the slope of the (longitudinal) pseudotime (or equivalently the 
longitudinal profile functions) is present at all values of the coupling and, indeed, was also visible 
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in the isotropic results (see Fig. 2 of Ref. |^]). This robustness suggests that this behaviour is not 
just an artefact of the various trial actions but, rather, that the variational calculation is attempting 
to mimic behaviour present in the exact theory. In other words, we suspect that if one were able 
to calculate the average mean square displacement in Eq. ( 2.24 ) with the exact weight function e*'^ 
rather than one would find equivalent behaviour. This expectation is supported by the fact that 
it was precisely this turnover in the profile function which was required in order to obtain, at least 
with the isotropic trial actions discussed in Sec. 4.3 of Ref. |Q, a non-vanishing total cross section at 
pion production threshold. 

Further evidence that anisotropy is quite important in this problem is shown in Fig. where the 
modified bare mass Mi is plotted as a function of the coupling constant a. Also shown are the results 
obtained with the increasingly sophisticated isotropic trial actions in Ref. (the results have been 
normalized by the ‘worst’ of these, namely Feynman’s ansatz for the retardation function in terms 
of an exponential, with two variational parameters). It is this quantity which serves as a ‘figure of 
merit’, because the variational principle guarantees that the true value of Mi is approached from 
below. We see that the improvement in this bound which is possible with the anisotropic trial action 
(as compared to that provided by Feynman’s parameterization) is essentially an order of magnitude 
larger than that obtained with the best isotropic action. Nevertheless, the values of Mi obtained 
in this work are at most only about 2 % above what one obtains with Feynman’s parameterization. 
This gives one the hope that the values of Mi shown in Table || and plotted in Fig. ^ may in fact be 
quite close to the exact ones (which, unfortunately, are not known at present but could be obtained 
through a direct numerical evaluation of the functional integral in Eq. (2.7)). We also stress, because 
the Ml’s obtained here are lower bounds, that this quantity serves as a useful yardstick to assess and 
compare the quality of the different nonperturbative calculations referred to in the Introduction. 

On the other hand, bearing in mind that in ordinary quantum mechanical applications of the 
variational principle even small gains in the ground state energy are usually associated with appreciable 
improvements in other observables - an experience which was also observed when Eeynman’s ansatz for 
the retardation function was replaced by the ‘improved’ or ‘extended’ (isotropic) parametrizations 
- one might also expect substantial improvement for form factors and amplitudes when the anisotropic 
rather than isotropic trial actions are employed. In particular, some unsatisfactory properties of the 
variational description of meson-nucleon scattering, such as the observed shift in the position of multi¬ 
meson thresholds away from s = (M -|- nm)"^, the violation of unitarity |^] and the unphysical tails in 
the scaling function Q are likely to be amended. 


Finally, we turn to the critical coupling of the (quenched) theory. As shown in Ref. Q, with 
Feynman’s trial action the variational equations cease to have real solutions above ac ~ 0.824. With 
the best isotropic action this reduces slightly to ac ~ 0.815. In the present work, with the most 
general anisotropic trial action, we find that this coupling reduces significantly to Uc ~ 0.709. It 
is remarkable that if one fixes A to the perturbative value 1 (i.e. if one removes it as a variational 
parameter), the anisotropy in the profile functions and pseudotimes is now sufficient to produce a 
critical coupling, whereas this didn’t happen in the isotropic case. Obviously the occurrence and the 
magnitude of the critical coupling are directly linked to the flexibility offered by the trial action and 
if it would be possible to go beyond quadratic trial actions one may expect that the critical coupling 
would decrease even further. It is tempting to speculate that the true ac (irrespective of whether this 
is 0 or positive) is being approached from above, although we know of no rigorous proof of this. The 
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Figure 3: The finite mass Mi as a function of a, normalized by the result obtained with the Feynman 
parameterization. Mi is a strict lower bound on the exact value. It is clearly seen that the anisotropic 
action leads to significant improvement as compared to the isotropic actions. These are labelled “Best” 
when no specific form for the retardation function has been assumed assuming and “Improved” for a 
parameterization which incorporates the correct behaviour at small times. 


critical coupling shows the same (almost) linear rise as a function of the pion mass which was observed 
in Ref. (see Fig. |^, with the results obtained with the anisotropic action lying about 15% below 
those obtained with the best isotropic action. On the other hand, the values of the critical coupling 
observed in the rainbow Dyson-Schwinger equation studies reported in Ref. |^] are somewhat larger 
than those labelled ’Best isotropic’ in Fig. 0, again with the same linear increase as a function of m. 


4 Summary and Conclusions 

We have demonstrated that the most general quadratic trial action in the polaron variational approach 
to the scalar Wick-Cutkosky model leads to a substantial improvement over the previous results which 
utilized only isotropic “profile functions”. By having a maximum principle for the intermediate mass 
Ml , this is clearly reflected in the values for this quantity displayed in Table |l|. Moreover, the critical 
value for the coupling constant a (i.e. the value of a above which no real solutions for the variational 
equations exist) obtained with the more general anisotropic trial action is lower than with the isotropic 
trial actions. This is consistent with the interpretation that the more general the trial action, the 
closer the critical coupling approaches its exact value ac = 0. Along the way we have introduced 
new, powerful, methods to evaluate various required functional averages, we have shown that the 
variational principle itself demands anisotropy when left enough freedom and discussed the physical 
picture for the different evolution of the dressed nucleon parallel or perpendicular to its momentum. 

Finally, it may be useful to discuss some distinctions and merits of the polaron worldline approach 
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Figure 4: The critical coupling as function of the pion mass, as obtained with the anisotropic action. 
For comparison, the results of Ref. Q for the isotropic trial actions are also shown (the so-called 
‘Feynman’ and ‘improved’ actions results were only calculated at the physical pion mass). The dotted 
curves have been drawn to guide the eye. 


compared to other non-perturbative methods. Despite its restriction to only quadratic trial actions 
this work has again shown that it is capable and flexible enough to cope with the many different 
scales occuring in a field theoretical problem. In addition, we would like to stress its truly variational 
aspect which allows to distinguish in a quantitative manner between the different ansdtze; in contrast, 
the so-called “variational” perturbation theory |26] is more an optimization of perturbation theory 
with respect to some artificially introduced parameter and relies on the ad hoc “principle of minimal 
sensitivity”. Schwinger-Dyson equations are widely used non-perturbative methods which, however, 
require truncations of the infinite hierarchy of equations connecting different Green functions. This 
truncation introduces uncertainties which are hard to control or to quantify and frequently leads to 
spurious gauge dependence of results when applied to gauge theories. In comparison, the worldline 
variational approach has been shown to respect gauge covariance [ p0| ] and can be corrected systemat¬ 
ically by calculating higher-order cumulants. Of course, nothing is known about the convergence of 
these corrections to the exact result (except that it seems to be rapid in the polaron case |2^) but 
this feature is shared by any other systematic non-perturbative method. Note also that no derivative 
expansion is introduced in the variational approach as it usually is in the “Exact Renormalization 
Group” approach (where it leads to a dependence on the regulator in every order of that scheme). 
Neither has the number of constituents to be small as in numerical applications of “Discretized Light- 
Gone Quantization” or the Tamm-Dancoff method |^. Again the polaron provides an instructive 
example that suggests that such a restriction may be only meaningful for the weak-coupling case: the 
mean number of phonons surrounding the dressed electron grows like 0.217a^ at large couplings |^]. 
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Of course, there are also many restrictions of the worldline variational method in its present form; 
perhaps the biggest deficiency is its limitation to the quenched approximation and abelian gauge the¬ 
ories which prevents application to realistic strong-coupling theories like Quantum Chromodynamics. 
Nevertheless, we believe that this approach has enough virtues to make further applications and 
extensions worthwhile. 


Acknowledgements: One of us (AWS) would like to thank Martin Oettel, Will Detmold and Alex 
Kalloniatis for numerous illuminating discussions. 
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Appendix 


A The behaviour of the quenched Wick-Cutkosky model at large 
orders of perturbation theory 

It is well known that the large order behaviour of a field theory’s perturbative expansion is a useful 


tool for ascertaining its stability |^]. In this Appendix we briefly summarize these arguments and 
apply them to the quenched cubic scalar theory considered in this paper. 

Half a century ago Dyson argued, on physical grounds, that the perturbative expansion 
Jk 9^ of typical Green functions, such as the propagator for the field ‘h which we consider here, 
should be asymptotic. Indeed, one typically finds that the size of the coefficients fx does grow 
factorially, preventing simple-minded summation of the series. Nevertheless, if subsequent terms in 
the expansion oscillate in sign or change in phase (as they do, for example, in <I>^ theory |^^), Borel 
re-summation techniques allow one to restore meaning to the expansion. Only when all terms at large 
K contribute with the same sign is this not possible. In this case, the perturbative vacuum is not the 


true vacuum and tunnelling to the true ground state occurs |^|. The typical, exponentially small, 
imaginary parts which this entails can be obtained by first moving the coupling g slightly off the real 
axis, therefore allowing the Borel re-summation of the perturbative expansion, and only at the end 
analytically continuing back to Im g = 0. In this way the divergent nature of an expansion generates 
imaginary parts even if fx is real at all orders. A good example of this is provided by unquenched 
theory which was studied with the techniques outlined here in Refs. (Earlier papers discussing 
the divergent nature of unquenched <I>^ theory may be found in Ref. [0^].) 

Dyson’s physical argument involves the different nature of pair production for positive and negative 
couplings and clearly is not relevant in a quenched theory such as the one considered here - the vacuum 
for the quenched theory is stable by construction. This need not be the case for the particle excitations, 
however. Indeed, for the propagator of one can make a rough estimate of fx by simply counting 
the number of Feynman diagrams at each order, namely 


f2K 


2^Kl 


V2K^ 


i + o|l 


(A.l) 


The typical factorial growth of these coefficients is present in the quenched theory just like in the 
unquenched one. Moreover, each order contributes with the same sign. This is readily seen from the 
worldline expression for the d-dimensional Euclidean coordinate space propagator 0 


G 2 {x) = const J dT exp 


• exp 



i'x{T)=x 

lx{0)=0 
n,2 nT 


Vx{t) 


+ dtidt 2 <x{ti)\Dj^\x{t 2 )>\ (A.2) 


Here < x{ti)\D^\x{t 2 ) > is the propagator for and is positive. In short, unless the behaviour in 
the above estimate for fx is grossly wrong, the nucleon will unavoidably be unstable at any coupling, 
with an exponentially small imaginary component to its mass. 

The estimate of fx may be refined through the standard functional techniques discussed in 
Refs. 1^ (see also Edwards |^| who, by different means, showed that the perturbation expansion of 
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quenched theory is asymptotic regardless of whether il) is fermionic or bosonic). The method we 
largely follow is that of Itzykson, Parisi and Zuber who discussed the vertex function in quenched 
QED. In fact, the present application is much simpler not only because there is no spin and no gauge 
symmetry to worry about, but also because in the quenched scalar theory we need not worry about 
renormalons or, indeed, UV divergences if we work in a suitably small spacetime dimension (as pointed 
out in Ref. |]R|, the instability should not be a function of dimension). 

The quenched propagator 


G2{x) = 


6J{x) 5J{0) 


J 'Dip'D^ exp|— y d'^x [C{x) + J(x)<l>(x)] 


quenched 


(A.3) 


j=o 


(£ is given in Eq. (|1.1D) is determined by the functional average of the propagator of ‘h in a background 
field generated by tp, weighted by the free action associated with ip: 


G2ix) = J Vip < X 


1 


y = 0 > exp 


J d‘^x ip{x) ip{x) 




(A.4) 


— □ + Mg — 2g(p 

where we have assumed that the functional integral is appropriately normalized. The largest contri¬ 
bution to the ip integral will come from the region where the denominator in Eq. (A.4) is smallest, 
i.e. when 

-O + - 2gip{x)] 4>(x) ~ 0 . (A.5) 

This suggests that we consider the ‘eigenvalue’ equation 


-a-\-- 2gn[(p] ipix) 4>„(x) = 0 , 


(A.6) 


where the index n enumerates the eigenvalues gn and eigenfunctions which both have functional 
dependence on ip. Clearly gn is antisymmetric, i.e. 


9n[-^] = 

and will turn out to be nonzero, while the appropriate orthogonality condition for is 

J d'^x^iix) ^j{x) ip{x) = 6ij J d'^x^i^x) ip{x) = 5ijNi[ip] . 

We make use of these considerations reverting to the functional integral over 4> in Eq. (| 
evaluating it by expanding 4> in terms of the solutions i.e. 

4>(x) = ^Cn$n , 

n 

the functional integral over <I> being replaced by integrals over all coefficients c^. We obtain 

4>n(0)$n(x) 


G2{x) = VipJ2- 


2 iVnb] (ffnM -g) 

and so the coefficient of g^^ is given by 


exp 


n 

.,2K ■ 


y J d'^x ip(x) □ -k ip(x) 


(A.7) 

(A.8) 
and 

(A.9) 


(A.IO) 


f f d>„(0)4>n(3;) 

f2K = Vip nAT l } -TT 

J ^2Nn[(p]gnm 


d'^x(p{x) (^-a + ip{x) - 2Kloggn[<p] 


. (A.ll) 
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Terms odd in K vanish because of the antisymmetry of gn- 

For large K the expression for f 2 K may be estimated by steepest descent in the usual way and 
hence one seeks solutions leading to finite action in Eq. (|A.ll), of 


/ n\ 2K 5 r -I 

(-□ + m^]ip + ——r —gn[^] 
V / gnM dip 


= 0 


(A.12) 


ip=ifiO 


The functional derivative of gn[p\ may be evaluated as in Ref. 

^9n[p\ _ gnW\^l{x) 


and is given by 


5(p{x) 


Nn[p\ 


(A.13) 


The dependence on K may be extracted from the coupled set of differential equations ( |A.6|) and (|A.13D 
by defining 




= Vk<Po, Nn = J d'^x ip{x) , gn = 


y/K 


9n 


and so the solutions to 


+ ^nix) - 2gn[po]ipo{x)^n{x) = 0 (po{x) 


d^n[Po] 


(A.14) 


<^l{x) = 0 . 


(A.15) 

are independent of AT. It is not difficult to convince oneself that normalizable solutions do exist: if 
one tries = \/Nn 9 n/K po one obtains the usual equation for the unquenched theory which has 
analytic solutions if d = 6 and m = 0 |35]. For d 7^ 6, m 7^ 0 and/or general <I> Eq. ( A.15| ) is readily 
solved numerically. 

One therefore obtains the large K behaviour of fx to be essentially the same as that estimated 
from diagrammatic counting, i.e. 


hK = 




Nn[po] 9n[po] {eglipol) 


K 




(A.16) 


Here detAf is the determinant resulting from the quadratic fluctuations around po- Care must 
be taken in its evaluation as the extraction of the zero modes associated with the translational and 
dilatation symmetries of Eq. (|A.13 ) introduce extra powers of we refer the reader to the extensive 
literature on this aspect |31]. The central point, however, is that the factor already seen in 
Eq. (A.l) remains. 


B Calculation of Averages 

Here we describe an efficient method to calculate the various averages needed in the application of 
the variational principle. It is based on functional integration over velocities, instead of coordinates, 
and a ‘Hilbert-space’ notation for the various proper-time dependent quantities which allows an easy 
evaluation of the T —> 00 - limit required for the on-mass-shell case. A preliminary account has been 
given in Ref. We will keep the space-time dimension d arbitrary in these appendices (thereby 

allowing utilization of the results for both the 3-dimensional polaron problem as well as dimensionally 
regularized theories), while in the main text d has been set to 4. 
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B.l Integration over velocities 

Integration over velocities offers some simplifications, in particular for the fermionic case |^]. In 
the bosonic case its main virtue is the absence of explicit boundary conditions and the simple form 
the most general quadratic trial action takes if expressed in velocity variables. Transition to these 
variables simply amounts to multiplying the discretized path integral 


N-l 


^ const j d‘^{Xb-Xa)e^P'^^’’ I , h-ta = iVAt, Xq = Xa, Xn = Xb 

(B.l) 


k=l 


by 


~ ^ n J d'^Vkd{xk-Xk-i-AtVk) . (B.2) 


After performing the x^-integrations (A; = 1,... A — 1) one obtains. 


^ pt 

Xk = Xa + At''^Vj , i.e. x{t) =Xa+ dt'v{t') . 

3=1 


(B.3) 


There is one 5-function left which expresses the constraint that full integration over the velocity has 
to give the final position 

x{tb) = Xb = Xa+ [ dt' v{t') (B.4) 


but this is removed by the integration over Xb — Xa- Thus 


Dx = const. / Dv exp ( iS x(t) = Xa+ dt' v{t') 


(B.5) 


The relation ( [B.3f ) singles out the initial point Xa where the particle starts. Using Eq. (B^) we may 
also use the final position Xb as reference point: x(t) = Xb — //*' dt' v{t') . This can be combined with 
Eq. (B.2) into the more symmetrical form 


x{t) = 


Xb + Xa 


+ 


I 


2 2 

where sgn(x) = 20(x) — I is the sign-function. 


'ta 


dt' sgn(t — t') v{t') 


(B.6) 


B.2 Hilbert space formulation 

It is useful to introduce a kind of Hilbert space notation and to write Eq. ( |B.6|) as 

x{t) = {t\x) = + i(t|g|„) ^ + i ^ *'(«|S|«')(*'K>) 

= + dt' S{t-t')v(t'). (B.7) 

We further define 

(t|o) = 1 (B.8) 


21 











so that the relation between coordinate and velocity can be written representation-free as 

\n = l») + ^-SIO (B.9) 

Here the sign operator S is defined to have the matrix elements 

(t|5'|t^) = sgn{t — t') (B.IO) 

and the inner product is defined by integration over the proper time from ta to t^. Note that 


o) = S\t = ta) = -S'|f = 4) 


(B.ll) 


and 


{t\x) = {x\t)* = ^ ^ ^ dt'v{t'){-)sgn{t-t') (B.12) 

since everything is real. Therefore 5^ = —S is an antihermitean operator. Since dsign{t — t')/dt = 
25{t — t') we have, of course, x{t) = v{t). Due to translational invariance we are also free to choose 
the initial and final points: Xa = 0, x;, = x and Xq = —x/2,Xfe = x/2 are most convenient. 

Let us first consider the previous (isotropic) trial actions in this new Hilbert space notation: the form 


(2.4) can now be written compactly as 


St = - 

y (x|x) + ZKq (a; l/diag “ /| x) 

(B.13) 

where 

fto+T 


(H l/diagi 4) 

:= 5{ti - t2) / dt' f{ti - t') . 

(B.14) 

J to 


Using the relation (|B.9D this becomes 


(B.15) 

where 



A = 


(B.16) 


is the “profile function”. The same expression for the velocity trial action is obtained from Eq. ( |2.5|) , 
the only difference being that here the relation between profile function A and retardation function g 
is 


A = l + 2iKog. (B.17) 

Actually, Eq. ( B.15 ) is also the most general quadratic action associated with Eqs. ( |2.10 , 2.11|) if we 
understand the profile function as a Lorentz tensor 


'S’tH = 


fLlA I 


Vu) , 


(B.18) 


i.e. in general the matrix A also has Lorentz indices. It is easy to derive its form from Eqs. (T mm 
but we may equally well decompose A from the available structures (the metric tensor and the tensor 
formed from the external momentum p), forgetting about the original retardation functions. For later 
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purposes it is most convenient to decompose (and other Lorentz tensors) into components parallel 
and perpendicular to p 


A^^it - t') = AUt - 0 ^ + Ar{t - t>) f,- - ^ 

pZ y pZ 


= AiPr + AtP^’' . 


(B.19) 


If no confusion arises we will not write the Lorentz indices explicitly in the following. The extended 
trial actions (2.9) then just add a linear term in v and may be written as 


Ko 


‘S'tN = - — {v\A\v) + {b\v) 


(B.20) 


where 

|5) = Ap|o). (B.21) 

Eq. ( B.20 ), together with Eqs. ( B.19 , |B.2lD , is the most general linear + quadratic trial action 


consistent with translation and time translation invariance. By construction it reduces to the free 
action for A = \ = 1. 


B.3 Averages 

For evaluation of the various averages needed in Eq. 
dimensions 


we use the simple Gaussian integral in d 


J V^v exp I i 


^ 2ko ^ 


-y(u|Al|u) + (6|u) 
^ 6^ + ^zTr log A 


const. 




=: e 


JF 


(B.22) 


as master integral. Here {t\b) is an arbitrary function needed to evaluate (^i). A special case is 
b{t) = \p needed for the average of St,So . The trace implies summation over Lorentz indices as well 


as integration over the continous proper time. Eq. (2.8) now becomes 

g™{p, T) = exp [ iPt - iPo + i (Sq - St + Si 


(B.23) 

where T), Fq refer to the case of trial and free (i.e. A = A = 1) action, respectively. Prom the master 


integral (|B.22 ) we obtain for the different averages 
a) Ft - Fq: 


Ft — Fn = 


1 

2ko 


b]-{bo\bo) 


+ —i Tr log A 


(B.24) 


where |5o) = p|o). 

b) {St): 

{St) = “^log J ex.p{irSt) 


r=l 


d 

= —Ft[rA,rb] 


r=l 


1 

2kq 


6 )+-iTr(l) (B.25) 


c) (^o): Since 6A{ti — t 2 )/bA{d) = 6{ti — t 2 ) we have 


(^o) = 


6 /■ , , / X <5 

+ dt bo{t) 


6A{0) 


6b{t) 


Ft = - b 

2ko 


A 2 


b]+— [b 

Ko 


bo]+pFrj. (B.26) 
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(B.27) 


d) (Si): According to Eq. ( p^) we have to work out 

exp { -ik ■ [x{ti) - x{t 2 )] } 

The above average may be evaluated again by means of the master integral (B.22|) with 


b{t) 


or representation-free 


We thus obtain 


bi{t) = Xp- - [sgn(ti -t) - sgn(t 2 - t) ] 


\bi) = \b) + -S[\h)-\t2) 


(B.28) 

(B.29) 


( exp{-z/c • [x(ti) - x(t 2 )] } ) = exp<j — 


hi 


h\-[b 


with 


Ul = 


2 

av f 

O2 = 7 

^ 4 L 


= exp 


1 S A*!/ 


2 ko . 


2 ai ■ k + k ■ a2 ■ k 




s 


h] - 


t2) 




=: E{k,ai,a2){B.30) 

(B.31) 

(B.32) 


In d = 4 — 2e dimensions one has to substitute where u is an arbitrary mass 

parameter, in order to keep a = g ^/dimensionless in any dimensions. We thus have 


1 


(Si) =-^ / dtidt 2 / 7 —w T 9 - 9 - 

^ ^ 2kI Jto J {27r)d k^ - + iO 


exp 


{ 2^ [ 2ai • /c -I- fc • 02 • A: ] I . 

(B.33) 


The /c-integration can be performed by exponentiating the meson propagator 

1 1 


/j.2 _ ^2 _|_ 2iKo 


du exp 


. 2 ko ^ 


^ I u' 


(B.34) 


Taking care of our -|- 


-metric in Minkowski space we obtain 


d-i 


is.) - 


2 kq 2iKo 

) 

du' 


{ 271 ^ 


'^0 


1 


0 deti/2(^Mu') 


exp 


2hq - 


m^u' + ai ■ C • oi 




with 


= u'g^^'' + of = {u' + + {u' + al)Pf . 


The decomposition into the orthogonal projectors of Eq. (|B.19 ) immediately gives 
1 1 


W')' 


pr + 


u' + 02 


L 


u' + O 2 


detc;; = {u + a^) ■ {u' + al) 


T\d-1 


(B.35) 

(B.36) 

(B.37) 
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Since ai oc p (see Eqs. ( B.3]j , [B.21 )) only the longitudinal part of 02 survives in the exponent. 
Substituting u = / {u' + a^) one then obtains 


(5i) = 


2i7ri/"^ 


IGvr^ACo I kq 




1 


rl 


I dtidt2 I 

ho (4) Jo 

1 -‘il2+e 


1 + M-- 1 U 


• exp ■ 


i 

2kq 


du u 


9 r 1 — u at 

m Oo-1—hu 

u ak 


. (B.38) 


Putting everything together we have 


—i 


loggTip^T) = — 
ZKq L 


ibo\bo)-2{b 


iT f 2 
2ko V 


1 

A 

-2A 


60 + & 


(o|Al 


-l| 


(o|o) 


A2 


+ A' 


1 


log^+ J- 1 


+ i(Si} (B.39) 


;(o|^ 


- 21 


(o|o) 


where 


n(T) = 


Ko 


2i(c 


■Tr 


log ^ ^ - 1 


V(T} = 

(op) 

As indicated, all these quantities are in general T-dependent. 


•p + 20(r) + 21/(r)|(B.40) 

(B.41) 

(B.42) 


C Variational equations for finite T 

Here we present the variational equations for the most general case: an unspecified profile function 
a free linear term b^ in the trial action and finite proper time. The latter case is needed, for 
example, for off-shell Green functions. Suppressing Lorentz indices and defining 


\b) = Al|c) 


Eq. (|B.39[) simply becomes 


log^^b;^) = (^ol^o) - 2 (c| 6 o) + (c|c) 

2Ko l 


- -Tr 
2 


logA + — - 1 


+ i{Si) 


According to Eq. (B.31), in the averaged interaction only 


ai = 1 1^1) “ 1^2) 


(C.l) 


(C.2) 


(C.3) 


depends on |c) whereas the sole dependence on A resides in 02 (see Eq. ( B.32 )). Instead of varying 
with respect to b{t) or the parameters contained therein we may vary equivalently with respect to c(t). 


As it displays the ai, a 2 -dependence in the simplest way, the average ( B.33 ) is the most convenient 
expression for the present purposes (Eq. ( p.38| ) is also applicable if longitudinal and transverse parts 
are taken with respect to the vector or c^). We then obtain the variational equation for c{t) as 

ho 


ci^{t) = 6K(i) + = 6(((t) - / (itidt 2 [sgn(t - ti) - sgn(t - t 2 )] 


(27r)‘^ A ;2 — m? + zO 


-E{k,ai,a2) (C.4) 
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where E{k,ai,a 2 ) is defined in Eq. ([B.3C|). Similarly one obtains the variational equation for the 


profile function by varying Eq. ( C.2 ) with respect to This gives 


= g^'^d(t — t') -/ dtidt2 [sgn(t — ti) — sgn(t — t2) ] 

8 kq J to 


[ sgn(t' - ti) - sgn{t' -t 2 )] J 


d^k 




Eqs. ( |C.4| , |C.5| ) are a system of coupled nonlinear integral equations which also determine the optimal 
Lorentz structure of the linear term and the profile function. Despite the complexity of these equations 
their very structure allows some general observations: first, it is seen that the profile function is 
symmetric A(t',t) = A{t,t'). Second, we note that both c and A take their free values at the proper 
time boundaries Iq, to+T, since then sgn(t — ti) — sgn(t —12) = 0. It is also seen that the inhomogenous 
term bglt) is essential for the Lorentz structure because it is the only available 4-vector. If it is absent 
(for example, in the case of the massless on-shell propagator) then = 0 and A^‘' oc . This is 
because the ^-integral in Eq. (C^) would be proportional to a((, i.e. and the /c-integral in Eq. 
(|C.5| ) proportional to a^ai and g^’^. Therefore, no anisotropy is generated in this case. By the same 
argument, a constant inhomogenous term 6 q (for example, in the case of the massive on- or off-shell 
propagator where 6 q = p^) naturally generates c^(t) = A(t)6Q and an anisotropic profile function like 
in Eq. ( B.19| ) where the Lorentz decomposition is with respect to the preferred vector bo. In the 
on-shell limit \{t) —> A and A{t',t) —> A{t — t') (except at the boundaries) and in d = 4 Euclidean 
dimensions we recover precisely the anisotropic variational equations (|2.26|) - (2.28). 


D Large-T limit 

To obtain the position of the pole in the nucleon propagator we have to work out the limit T —> oo 
for the various averages calculated above. This is done most conveniently by choosing the symmetric 
proper time interval, i.e. to = —Tj^. 


D.l Leading order 

For large T all quantities 0{t — t') which only depend on the time difference may be diagonalized in 
the space of normalized functions 


{t\E) := 


I 


AEt 




(D.l) 


VIZ. 


{E\0\E') 


r+T/2 

/ dtdt' 0{t — t') ex.p{—iEt + iE't') 
J-T/2 


T—>O0 


27r J-T/2 

1 r+’^ f E + E' \ rHT-W\)/‘2 

— daO{a)exp(-i —-— a) j 

ZTT J-T \ 2 J J_(T-|o-|)/2 

r+oo 

6{E - E') / da 0{a) =: 5{E - E') d{E) . 


dE (D.2) 


(D.3) 


For notational simplicity we suppress the “tilde”-sign over the Fourier-transformed quantities in the 
following, i.e. write A{E) for A{E) etc. Using 

(E|o) = V^5{E) (D.4) 
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S{E) = -2iV^, (D.5) 

where V denotes the principal value, one can then easily evaluate the required large-T limits. For 
example 


1 y' 

A 


dtdt' ^ r°°dE ( 


J-T/2 

1 r+oo 

27r J —rvi 


27r J—oo 

pT riT-\a\)/2 

dE da dE e*®'" 


l-T 


l-{T-\a\)/2 


1 T- 

A{E)) 


/I \ 

^°°r(^TVTl .(D.6) 


U(o); 


The last equation can be obtained more easily by using Eq. & which gives 27r<5(0)/^(0) and 
replacing 27r(5(0) by T. This is similar to calculating a scattering cross section from the square of 
a transition matrix which contains an energy-conserving <5-function. The same rule applies if one 
evaluates (o|o) = T in the energy representation or {E\E') = 5{E — E') —> T/(27r) for E —> E'. 

Using the traces over Lorentz indices 


trP l = 1 , trPr = = d-1 


p. 


(D.7) 


one then obtains for the ‘kinetic term’ (B.41) 


0 = hm^tr/ dE{E\E) 
T >oo ZitJ J —oo 

= ^n[AL] + ^n[AT] 


r+oo 


log A{E) + 


A{E) 


- 1 


Al{E) 


— 1 -|- (d — 1) 


log At {E) + 


At{E) 


- 1 


(D.8) 


where we have used that A(E) is even. This property follows from the connection ( [B.16| ) in P-space 

(D.9) 


A(E)6(E-E') = S(E-E') + ^Ko^(Elf-U,,lE')^ 


or 


ME) = l + ‘^^ J_ da f{a) 


^iEa _ ^ 


= 1 - 


SIkq 

~w 


da f{a) sin" 


(f) 


(D.IO) 


Since the retardation function f{a = p — 12) must be even, it follows that A{E) = A{—E). The same 
conclusion is reached if the equivalent form (|B.17) is used: 


roo 

A{E) = 1 + Aikq / da g{a) cos{Ea) . 
Jo 


(D.ll) 


Note that Eqs. ( p.lO , D.ll) are consistent with the relation f{a) = g{a). 

Next we evaluate the coefficients in the exponent of the averaged interaction term for large T and 


we obtain by using Eqs. (|B.3l|, |B.32| , p.5[) 

T^oc 1 


a!i{a) 


2 7-00 
= -a\p^ 


dE 6{E) A m 


1 ^-2i 


AiE)) 


V- 


( 1 


Am 


E 

p'' = —aXp 




(D.12) 
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»2 {<^) -I -iEV—g^ 

s\n^Ea/2 ( 1 y _ 


Stt J- 
4 
vr Jo 


oo 

oo „;„2 


giSil _ gji?t2 










— El{^) Pl^ + Pt^ 


(D.13) 


Here we have used the decomposition oil/A into orthogonal projectors Pi,Px and the abbreviation 
( 2.16|) . Finally, in the limit T —> oo the double integral in < S*! > simplifies to 


rT/2 


l-T/2 


dtidt2 


da 


r+(T-\a\)/2 


dS ... 


T—^oo 


P + OO 


’-T J-(T-\a\)/2 


da (T — |(t|) ... = 2 / da {T — a)... 
Jo 

(D.14) 

since the integrand does not depend on E = (H + t2)/2 and is even in a. The potential term ( |B.42D 
therefore becomes 


V 


T—^OO 


(2i'Kv‘^\ 


da 


87r2 \ Ko J Jo [glia)]^ 


■ exp 


2ko 


du u 


m^gl{a) 


1 + 


Vr(c^) 

1 — rt X^p 


- 1 U 


-3/2+e 


+ 




U 




Thus Eq. ( B.4C| ) has the following large-T limit 

i 

2ko - 


log T) - — Ip"^ {I- Xf + 2n + 2v] T + O M 

2 kq L J V / 


(D.15) 


(D.16) 


D.2 Subasymptotic terms 

The ©(T'^j-terms are needed for evaluation of the residue of the propagator. They can be obtained 
from Eq. (D.2) which for hnite T reads 

1 /■+'^ 


{E\0\E') = ^ daO{a) 


E + E' \ sm[{E - E'){T -\a\)/2] 


exp [ —-—a 


E- E' 


(D.17) 


Since with y = {T — |ct|)/2 


'■+°° sin(xy) 


= /(O) + /■+“ MuM 

X J-oo X J-oo X 

7r/(0)-2® 


(D.18) 


one obtains 


{E\0\E') = 5{E-E')-^5'{E-E') + ... da 0{a) exp (^-i^P^a'^ . (D.19) 

'-V-' 

=:OTiiE+E')/2) 

As expected, in next-to-leading order the operator O is no longer diagonal in E-space. More generally, 
any function of the operator O has matrix elements 


(E|E(0)|EO = 5{E-E')E{OTiE))-^6'{E-E')OT(^ 


E + E' 


E' (Ot 


E + E' 


))+... (D.20) 
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This allows one to evaluate the corrections for large T where we assume that the integration limits in 
Ot can be extended to ±oo with impunity (it should be remembered that the retardation functions 
and — t') decay exponentially for large time differences). However, due to the (^'-function these 

corrections involve derivatives of even functions and therefore lead to a vanishing contribution. For 


example, the kinetic term (B.41) 


H(T) = tr 

^ ^ 2iT 


C + OQ 


r+T/2 


I dEclE' I dt {E 

1-00 J-T/2 


log -4 + — - 1 


E'\ {E'\t) {t\E) 


+ 00 


dEdE' E 


■^tr / 

2iT J-oo 


does not receive a subasymptotic correction: 

n(T) - n dy 1- 


log H + ^ - 1 


1 


E' 


1 sin[(F; - E')T/2] 


TT 


E-E' 


r^dx6'(x) = 0 

A{y) J J-oo ^ ’ xT/2 


(D.21) 


(D.22) 


Similarly the corrections to other quantities either lead to odd integrands or to terms involving H'(0). 
Therefore all these corrections vanish and the only contribution in next-to-leading order comes from 
Eq. (pT4|) 


ViT) 


T^oo 


1 

V -\ - — 

T 87r2 


2mv^ 

Kq 


da ■ 


a 




l-e 


duu 


1 + 


AtT(o') 


-1 -3/2-re 


- 1 u 


• exp 


2kq 


m y,^{a) 


I - u X^p^a"^ 


+ 


u 


ni^) 


(D.23) 


E Virial theorem 

Here we prove that the ‘kinetic term’ 

H = (u • 


dE 


log A{E) + 


A{E) 


- 1 


(E.l) 


may be expressed in terms of the (unspecified) ‘potential’ V[fjp{a)]. This relation only holds if the 
variational equation 

5(H + 1/) = 0 ^A{E) = 1--H2(E)-|^ (E.2) 


is fulfilled. The value of the constant cu will turn out to be irrelevant and for simplicity we consider 
here only the isotropic case, the generalization to the anisotropic one being straightforward. 


We first note from the definition (2.22) of the pseudotime that 

5y?{a) ‘isw?{Ea/2) 


6A{E) 

and therefore by multiplication with A(E) 


TT E‘^A^{E) 




(E.3) 


(E.4) 
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Next we split up O = fli + ^2 into two parts and evaluate them separately. We write 


rii = uj ■ / dE 
Jo 


1 


- 1 


(E.5) 


[A{E) 

and insert the variational solution ([E.2|) divided by A{E) into that expression. In this way one obtains 


poo POO POO 


5V 5ix‘^{a) 


6A{E) 

r ^ 2/ ^ 

/ (^ K 2( \ 


0 < 5^2 


(E.6) 


where Eq. (E.4) has been used in the last line. After an integration by parts the logarithmic term 
becomes 

n, = -..faEE^. (E,7) 

Here we have assumed for simplicity that the boundary terms do not give a contribution but it can 
be shown that the final outcome is the same even if this is not the case. H 2 can be transformed by 
using the variational equation differentiated with respect to the variable E 


d fsm\Ea/2y 

LOTT Jo oE y 


E2 


(E.8) 


Converting the derivative with respect to E into one with respect to a one obtains 


EA'{E) 


da 




, (9 / 1 4 sm'^{Ea/2) 
UJ Jo 6fi‘^{a) da icr^yr E‘^A{E) 


and therefore 


= - 


da 


Sy 3_^ ( 

5iJ?{a) da I a"^ 


Combining both expressions we have 


(5E 


Jo Sfi^ia)'^ da 


■ cr 


a ff^Ha)' 


a 


roo XT/ 

da „ r., . E^{a)—a 


5^j?{a) 


da 


(E.9) 


(E.IO) 


(E.ll) 


It is instrnctive to compare this resnlt with the usual virial theorem in nonrelativistic quantum 
mechanics which states that for one particle moving in a potential V{x) the expectation value of 
the kinetic energy is given by {xV'{x)) /2. Since in lowest order fjJ{a) = a one also sees that the 
variational kinetic energy is proportional to (conpling constant)^ for small coupling. Einally, it should 
be noted that the proportionality constant u) does not appear in Eq. ( [E.ll[) . Therefore the result can 
be taken over immediately to the anisotropic trial action and gives Eq. O- 
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